Detailed Seed Cone Morpho-Anatomy Provides New Insights into Seed Cone Origin and Evolution of Podocarpaceae; Podocarpoid and Dacrydioid Clades.

The study of reproductive morphology and trait evolution provides a vital insight to understand the evolutionary history of plants. The conifer family Podocarpaceae has a remarkable diversity of seed cones, with distinct morphology among the genera and with conifers in general. However, we lack a good understanding of the seed cone morpho-anatomy and trait evolution of Podocarpaceae. We investigated detailed seed cone morpho-anatomy using staining and sectioning techniques to clarify the anatomical, morphological diversity and evolution of functional traits. The presence of a fleshy receptaculum is a characteristic feature of both clades. However, species of Retrophyllum, Afrocarpus and some species of Nageia and Podocarpus form a fleshy sarcotesta-like seed coat, lacking a fleshy receptaculum. The ancestral state reconstructions show a shift between and sometimes within the genus. Although both clades demonstrate fleshiness as an ancestral trait, the shift in fleshy structures provides evidence for complex multiple evolutions of fleshy morphologies. These seed cone traits (e.g., fleshiness and size), along with the broad, flattened and well-adapted (leaf dimorphism) foliage in both clades, are largely congruent with efficient light harvesting and bird dispersal. These traits make these two clades well adapted to their environment, when growing in communities including tall and broad-leaved angiosperms (closed-canopy angiosperm forests), compared to other podocarps, making them more successful in achieving a wider distribution and species richness.

Diversity, Distribution, Systematics and Conservation Status of Podocarpaceae.

Among conifer families, Podocarpaceae is the second largest, with amazing diversity and functional traits, and it is the dominant Southern Hemisphere conifer family. However, comprehensive studies on diversity, distribution, systematic and ecophysiological aspects of the Podocarpaceae are sparse. We aim to outline and evaluate the current and past diversity, distribution, systematics, ecophysiological adaptations, endemism, and conservation status of podocarps. We analyzed data on the diversity and distribution of living and extinct macrofossil taxa and combined it with genetic data to reconstruct an updated phylogeny and understand historical biogeography. Podocarpaceae today contains 20 genera and approximately 219 taxa (201 species, 2 subspecies, 14 varieties and 2 hybrids) placed in three clades, plus a paraphyletic group/grade of four distinct genera. Macrofossil records show the presence of more than 100 podocarp taxa globally, dominantly from the Eocene-Miocene. Australasia (New Caledonia, Tasmania, New Zealand, and Malesia) is the hotspot of living podocarps diversity. Podocarps also show remarkable adaptations from broad to scale leaves, fleshy seed cones, animal dispersal, shrubs to large trees, from lowland to alpine regions and rheophyte to a parasite (including the only parasitic gymnosperm-Parasitaxus) and a complex pattern of seed and leaf functional trait evolution.

Population decline in a Pleistocene refugium: Stepwise, drought-related dieback of a South Australian eucalypt.

Refugia can facilitate the persistence of species under long-term environmental change, but it is not clear if Pleistocene refugia will remain functional as anthropogenic climate change progresses. Dieback in populations restricted to refugia therefore raises concerns about their long-term persistence. Using repeat field surveys, we investigate dieback in an isolated population of Eucalyptus macrorhyncha during two droughts and discuss prospects for its continued persistence in a Pleistocene refugium. We first confirm that the Clare Valley in South Australia has constituted a long-term refugium for the species, with the population being genetically highly distinct from other conspecific populations. However, the population lost >40 % of individuals and biomass through the droughts, with mortality being just below 20 % after the Millennium Drought (2000-2009) and almost 25 % after the Big Dry (2017-2019). The best predictors of mortality differed after each drought. While north-facing aspect of a sampling location was significant positive predictor after both droughts, biomass density and slope were significant negative predictors only after the Millennium Drought, and distance to the north-west corner of the population, which intercepts hot, dry winds, was a significant positive predictor after the Big Dry only. This suggests that more marginal sites with low biomass and sites located on flat plateaus were more vulnerable initially, but that heat-stress was an important driver of dieback during the Big Dry. Therefore, the causative drivers of dieback may change during population decline. Regeneration occurred predominantly on southern and eastern aspects, which would receive the least solar radiation. While this refugial population is experiencing severe decline, some gullies with lower solar radiation appear to support relatively healthy, regenerating stands of red stringybark, providing hope for persistence in small pockets. Monitoring and managing these pockets during future droughts will be essential to ensure the persistence of this isolated and genetically unique population.

Detailed seed cone morpho-anatomy of the Prumnopityoid clade: an insight into the origin and evolution of Podocarpaceae seed cones.

BACKGROUND AND AIMS: Seed cone traits are significant for understanding the evolutionary history of conifers. Podocarpaceae has fleshy cones with a distinct morphology compared with other conifers. However, we have a poor understanding of the seed cone morphology of the Prumnopityoid clade and within Podocarpaceae. This study presents detailed seed cone morpho-anatomy and the evolution of fleshy structures traits in the Prumnopityoid clade. METHODS: We investigated the detailed seed cone morpho-anatomy of selected species from the nine genera using the histological method. The evolution of morpho-anatomical traits was assessed using ancestral state reconstruction methods. KEY RESULTS: The Prumnopityoid clade has evolved fleshy seed cones using different functional structures (e.g. aril, epimatium, bracts or receptaculum) and fleshiness is an ancestral trait in the clade. An epimatium is present in all genera except Phyllocladus, but with different structural morphologies (e.g. a fleshy asymmetrical cup-like epimatium or an epimatium that is fused with the integument, forming a fleshy sarcotesta-like seed coat). In all species with fleshy sarcotesta-like seed coats, the endotesta is hard and woody, forming a sclerotesta-like structure and the epimatium and exotesta are fused, forming a fleshy sarcotesta-like structure. CONCLUSIONS: This study highlights that the Prumnopityoid clade has an amazing diversity of structures and complex evolutionary patterns. Fleshiness is an ancestral trait of the clade and has been achieved via diverse evolutionary pathways and structures. This clade has four distinct seed cone types, i.e. drupe-like, receptaculate, arilloid and dacrydioid cones, based on morpho-anatomical structures and traits. The macrofossil record also demonstrates the presence of several structures and traits.

A targeted capture approach to generating reference sequence databases for chloroplast gene regions.

Metabarcoding has improved the way we understand plants within our environment, from their ecology and conservation to invasive species management. The notion of identifying plant taxa within environmental samples relies on the ability to match unknown sequences to known reference libraries. Without comprehensive reference databases, species can go undetected or be incorrectly assigned, leading to false-positive and false-negative detections. To improve our ability to generate reference sequence databases, we developed a targeted capture approach using the OZBaits_CP V1.0 set, designed to capture chloroplast gene regions across the entirety of flowering plant diversity. We focused on generating a reference database for coastal temperate plant species given the lack of reference sequences for these taxa. Our approach was successful across all specimens with a target gene recovery rate of 92%, which was achieved in a single assay (i.e., samples were pooled), thus making this approach much faster and more efficient than standard barcoding. Further testing of this database highlighted 80% of all samples could be discriminated to family level across all gene regions with some genes achieving greater resolution than others-which was also dependent on the taxon of interest. Thus, we demonstrate the importance of generating reference sequences across multiple chloroplast gene regions as no single loci are sufficient to discriminate across all plant groups. The targeted capture approach outlined in this study provides a way forward to achieve this.

Phylogenomics and continental biogeographic disjunctions: insight from the Australian starflowers (Calytrix).

PREMISE: Continental-scale disjunctions and associated drivers are core research interests in biogeographic studies. Here, we selected a species-rich Australian plant genus (Calytrix; Myrtaceae) as a case study to investigate these patterns. Species of this endemic Australian starflower genus have a disjunct distribution across the mesic fringes of the continent and are largely absent from the arid center. METHODS: We used high-throughput sequencing to generate unprecedented resolution and near complete species-level nuclear and plastid phylogenies for Calytrix. BioGeoBEARS and biogeographic stochastic mapping were used to infer ancestral areas, the relative contributions of vicariance and dispersal events, and directionality of dispersal. RESULTS: Present-day disjunctions in Calytrix are explained by a combination of scenarios: (1) retreat of multiple lineages from the continental center to the more mesic fringes as Australia became progressively more arid, with subsequent extinction in the center as well as (2) origination of ancestral lineages in southwestern Australia (SWA) for species-rich clades. The SWA biodiversity hotspot is a major diversification center and the most common source area of dispersals, with multiple lineages originating in SWA and subsequently spreading to the adjacent arid Eremaean region. CONCLUSIONS: Our results suggest that major extinction, as a result of cooling and drying of the Australian continent in the Eocene-Miocene, shaped the present-day biogeography of Calytrix. We hypothesize that this peripheral vicariance pattern, which is similar to the African Rand flora, may explain the disjunctions of many other Australian plant groups. Further studies with densely sampled phylogenies are required to test this hypothesis.

Historical biogeography of Pomaderris (Rhamnaceae): Continental vicariance in Australia and repeated independent dispersals to New Zealand.

AIM: Gondwanan biogeographic patterns include a combination of old vicariance events following the breakup of the supercontinent, and more recent long-distance dispersals across the southern landmasses. Floristic relationships between Australia and New Zealand have mostly been attributed to recent dispersal events rather than vicariance. We assessed the biogeographic history of Pomaderris (Rhamnaceae), which occurs in both Australia and New Zealand, by constructing a time-calibrated molecular phylogeny to infer (1) phylogenetic relationships and (2) the relative contributions of vicariance and dispersal events in the biogeographic history of the genus. LOCATION: Australia and New Zealand. METHODS: Using hybrid capture and high throughput sequencing, we generated nuclear and plastid data sets to estimate phylogenetic relationships and fossil calibrated divergence time estimates for Pomaderris. BioGeoBEARS and biogeographical stochastic mapping (BSM) were used to assess the ancestral area of the genus and the relative contributions of vicariance vs dispersal, and the directionality of dispersal events. RESULTS: Our analyses indicate that Pomaderris originated in the Oligocene and had a widespread Australian distribution. Vicariance of western and eastern Australian clades coincides with the uplift of the Nullarbor Plain c. 14 Ma, followed by subsequent in-situ and within-biome diversification with little exchange across regions. A rapid radiation of southeastern Australian taxa beginning c. 10 Ma was the source for at least six independent long-distance dispersal events to New Zealand during the Pliocene-Pleistocene. MAIN CONCLUSIONS: Our study demonstrates the importance of dispersal in explaining not only the current cross-Tasman distributions of Pomaderris, but for the New Zealand flora more broadly. The pattern of multiple independent long-distance dispersal events for Pomaderris, without significant radiation within New Zealand, is congruent with other lowland plant groups, suggesting that this biome has a different evolutionary history compared with the younger alpine flora of New Zealand, which exhibits extensive radiations often following single long distance dispersal events.

Extinction pulse at Eocene-Oligocene boundary drives diversification dynamics of two Australian temperate floras.

The diversification dynamics of the Australian temperate flora remains poorly understood. Here, we investigate whether differences in plant richness in the southwest Australian (SWA) biodiversity hotspot and southeast Australian (SEA) regions of the Australian continent can be attributed to higher net diversification, more time for species accumulation, or both. We assembled dated molecular phylogenies for the 21 most species-rich flowering plant families found across mesic temperate Australia, encompassing both SWA and SEA regions, and applied a series of diversification models to investigate responses across different groups and timescales. We show that the high richness in SWA can be attributed to a higher net rate of lineage diversification and more time for species accumulation. Different pulses of diversification were retrieved in each region. A decrease in diversification rate across major flowering plant lineages at the Eocene-Oligocene boundary (ca 34 Ma) was witnessed in SEA but not in SWA. Our study demonstrates the importance of historical diversification pulses and differential responses to global events as drivers of present-day diversity. More broadly, we show that diversity within the SWA biodiversity hotspot is not only the result of recent radiations, but also reflects older events over the history of this planet.

Identifying fossil Myrtaceae leaves: the first described fossils of Syzygium from Australia.

PREMISE OF THE STUDY: Although leaves of Myrtaceae are easily identified to family level, very few studies have convincingly identified fossil Myrtaceae leaves to living genera. We used a broadly comparative approach with a large data set of extant taxa to confidently assign the mummified remains of myrtaceous leaves from early Miocene sediments at Kiandra (New South Wales, Australia) to a living genus. METHODS: Fossils were identified using a nearest living relative approach, against a database of 232 extant broadleaf rainforest species of Myrtaceae. Leaf cuticles were prepared from 106 species, sourced from herbarium specimens as well as some living individuals, and a further 127 records were assembled from the literature. A set of simple but phylogenetically informative cuticular characters were observed, described, and recorded under both scanning electron microscopy and standard light microscopy. KEY RESULTS: A new fossil species of Syzygium Gaertn. is described from mummified remains found in early Miocene (21.5-21.7 Ma) sediments. The fossil taxon is here named Syzygium christophelii sp. nov., in honor of the late Australian paleobotanist David Christophel. CONCLUSIONS: These fossils represent some of the most confidently described Myrtaceae leaf fossils published to date and are the first and oldest described fossil record of Syzygium from Australia. While several fossil parataxa have been illustrated from New Zealand, and several fossil species of Syzygium have previously been proposed in the literature, many of these fossils lack characters for a confident diagnosis.

Targeted capture to assess neutral genomic variation in the narrow-leaf hopbush across a continental biodiversity refugium.

The Adelaide geosyncline, a mountainous region in central southern Australia, is purported to be an important continental refugium for Mediterranean and semi-arid Australian biota, yet few population genetic studies have been conducted to test this theory. Here, we focus on a plant species distributed widely throughout the region, the narrow-leaf hopbush, Dodonaea viscosa ssp. angustissima, and examine its genetic diversity and population structure. We used a hybrid-capture target enrichment technique to selectively sequence over 700 genes from 89 individuals across 17 sampling locations. We compared 815 single nucleotide polymorphisms among individuals and populations to investigate population genetic structure. Three distinct genetic clusters were identified; a Flinders/Gammon ranges cluster, an Eastern cluster, and a Kangaroo Island cluster. Higher genetic diversity was identified in the Flinders/Gammon Ranges cluster, indicating that this area is likely to have acted as a refugium during past climate oscillations. We discuss these findings and consider the historical range dynamics of these populations. We also provide methodological considerations for population genomics studies that aim to use novel genomic approaches (such as target capture methods) on non-model systems. The application of our findings to restoration of this species across the region are also considered.

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species.

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow-leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid-capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single-nucleotide polymorphisms (SNPs) for FST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.

Identifying Centres of Plant Biodiversity in South Australia.

We aimed to identify regional centres of plant biodiversity in South Australia, a sub-continental land area of 983,482 km2, by mapping a suite of metrics. Broad-brush conservation issues associated with the centres were mapped, specifically climate sensitivity, exposure to habitat fragmentation, introduced species and altered fire regimes. We compiled 727,417 plant species records from plot-based field surveys and herbarium records and mapped the following: species richness (all species; South Australian endemics; conservation-dependent species; introduced species); georeferenced weighted endemism, phylogenetic diversity, georeferenced phylogenetic endemism; and measures of beta diversity at local and state-wide scales. Associated conservation issues mapped were: climate sensitivity measured via ordination and non-linear modelling; habitat fragmentation represented by the proportion of remnant vegetation within a moving window; fire prone landscapes assessed using fire history records; invasive species assessed through diversity metrics, species distribution and literature. Compared to plots, herbarium data had higher spatial and taxonomic coverage but records were more biased towards major transport corridors. Beta diversity was influenced by sampling intensity and scale of comparison. We identified six centres of high plant biodiversity for South Australia: Western Kangaroo Island; Southern Mount Lofty Ranges; Anangu Pitjantjatjara Yankunytjatjara lands; Southern Flinders Ranges; Southern Eyre Peninsula; Lower South East. Species composition in the arid-mediterranean ecotone was the most climate sensitive. Fragmentation mapping highlighted the dichotomy between extensive land-use and high remnancy in the north and intensive land-use and low remnancy in the south. Invasive species were most species rich in agricultural areas close to population centres. Fire mapping revealed large variation in frequency across the state. Biodiversity scores were not always congruent between metrics or datasets, notably for categorical endemism to South Australia versus georeferenced weighted endemism, justifying diverse approaches and cautious interpretation. The study could be extended to high resolution assessments of biodiversity centres and cost:benefit analysis for interventions.

Transcriptome sequencing, annotation and polymorphism detection in the hop bush, Dodonaea viscosa.

BACKGROUND: The hop bush, Dodonaea viscosa, is a trans-oceanic species distributed oversix continents. It evolved in Australia where it is found over a wide range of habitat types and is an ecologically important species. Limited genomic resources are currently available for this species, thus our understanding of its evolutionary history and ecological adaptation is restricted. Here, we present a comprehensive transcriptome dataset for future genomic studies into this species. METHODS: We performed Illumina sequencing of cDNA prepared from leaf tissue collected from seven populations of D. viscosa ssp. angustissima and spatulata distributed along an environmental gradient in South Australia. Sequenced reads were assembled to provide a transcriptome resource. Contiguous sequences (contigs) were annotated using BLAST searches against the NCBI non-redundant database and gene ontology definitions were assigned. Single nucleotide polymorphisms were detected for the establishment of a genetic marker set. A comparison between the two subspecies was also carried out. RESULTS: Illumina sequencing returned 268,672,818 sequence reads, which were de novoassembled into 105,125 contigs. Contigs with significant BLAST alignments (E value < 1e(-5))numbered at 44,191, with 38,311 of these having their most significant hits to sequences from land plant species. Gene Ontology terms were assigned to 28,440 contigs and KEGG analysis identified 146 pathways that the gene products from 5,070 contigs are potentially involved in. The subspecies comparison identified 8,494 fixed SNP differences across 3,979 contiguous sequences, indicating a level of genetic differentiation between them. Across all samples, 248,235 SNPs were detected. CONCLUSIONS: We have established a significant genomic data resource for D. viscosa,providing a comprehensive transcriptomic reference. Genetic differences among morphologically distinct subspecies were found. A wide range of putative gene regions were identified along with a large set of variable SNP markers, providing a basis for studies into the evolution and ecological adaptation of D. viscosa.

Breakdown of phylogenetic signal: a survey of microsatellite densities in 454 shotgun sequences from 154 non model eukaryote species.

Microsatellites are ubiquitous in Eukaryotic genomes. A more complete understanding of their origin and spread can be gained from a comparison of their distribution within a phylogenetic context. Although information for model species is accumulating rapidly, it is insufficient due to a lack of species depth, thus intragroup variation is necessarily ignored. As such, apparent differences between groups may be overinflated and generalizations cannot be inferred until an analysis of the variation that exists within groups has been conducted. In this study, we examined microsatellite coverage and motif patterns from 454 shotgun sequences of 154 Eukaryote species from eight distantly related phyla (Cnidaria, Arthropoda, Onychophora, Bryozoa, Mollusca, Echinodermata, Chordata and Streptophyta) to test if a consistent phylogenetic pattern emerges from the microsatellite composition of these species. It is clear from our results that data from model species provide incomplete information regarding the existing microsatellite variability within the Eukaryotes. A very strong heterogeneity of microsatellite composition was found within most phyla, classes and even orders. Autocorrelation analyses indicated that while microsatellite contents of species within clades more recent than 200 Mya tend to be similar, the autocorrelation breaks down and becomes negative or non-significant with increasing divergence time. Therefore, the age of the taxon seems to be a primary factor in degrading the phylogenetic pattern present among related groups. The most recent classes or orders of Chordates still retain the pattern of their common ancestor. However, within older groups, such as classes of Arthropods, the phylogenetic pattern has been scrambled by the long independent evolution of the lineages.

Leaf evolution in Southern Hemisphere conifers tracks the angiosperm ecological radiation.

The angiosperm radiation has been linked to sharp declines in gymnosperm diversity and the virtual elimination of conifers from the tropics. The conifer family Podocarpaceae stands as an exception with highest species diversity in wet equatorial forests. It has been hypothesized that efficient light harvesting by the highly flattened leaves of several podocarp genera facilitates persistence with canopy-forming angiosperms, and the angiosperm ecological radiation may have preferentially favoured the diversification of these lineages. To test these ideas, we develop a molecular phylogeny for Podocarpaceae using Bayesian-relaxed clock methods incorporating fossil time constraints. We find several independent origins of flattened foliage types, and that these lineages have diversified predominantly through the Cenozoic and therefore among canopy-forming angiosperms. The onset of sustained foliage flattening podocarp diversification is coincident with a declining diversification rate of scale/needle-leaved lineages and also with ecological and climatic transformations linked to angiosperm foliar evolution. We demonstrate that climatic range evolution is contingent on the underlying state for leaf morphology. Taken together, our findings imply that as angiosperms came to dominate most terrestrial ecosystems, competitive interactions at the foliar level have profoundly shaped podocarp geography and as a consequence, rates of lineage diversification.

Evolution of exceptional species richness among lineages of fleshy-fruited Myrtaceae.

BACKGROUND AND AIMS: The angiosperm family Myrtaceae comprises 17 tribes with more than half of the estimated 5500 species being referred to the fleshy-fruited and predominantly rainforest associated Syzygieae and Myrteae. Previous studies suggest that fleshy fruits have evolved separately in these lineages, whereas generally shifts in fruit morphology have been variously implicated in diversification rate shifts among angiosperms. A phylogenetic hypothesis and estimate divergence times for Myrtaceae is developed as a basis to explore the evidence for, and drivers of, elevated diversification rates among the fleshy-fruited tribes of Myrtaceae. METHODS: Bayesian phylogenetic analyses of plastid and nuclear DNA sequences were used to estimate intertribal relationships and lineage divergence times in Myrtaceae. Focusing on the fleshy-fruited tribes, a variety of statistical approaches were used to assess diversification rates and diversification rate shifts across the family. KEY RESULTS: Analyses of the sequence data provide a strongly supported phylogenetic hypothesis for Myrtaceae. Relative to previous studies, substantially younger ages for many of the clades are reported, and it is argued that the use of flexible calibrations to incorporate fossil data provides more realistic divergence estimates than the use of errorless point calibrations. It is found that Syzygieae and Myrteae have experienced elevated diversification rates relative to other lineages of Myrtaceae. Positive shifts in diversification rate have occurred separately in each lineage, associated with a shift from dry to fleshy fruit. CONCLUSIONS: Fleshy fruits have evolved independently in Syzygieae and Myrteae, and this is accompanied by exceptional diversification rate shifts in both instances, suggesting that the evolution of fleshy fruits is a key innovation for rainforest Myrtaceae. Noting the scale dependency of this hypothesis, more complex explanations may be required to explain diversification rate shifts occurring within the fleshy-fruited tribes, and the suggested phylogenetic hypothesis provides an appropriate framework for this undertaking.

Comparative study of the evolution of nuclear ribosomal spacers incorporating secondary structure analyzes within Dodonaeoideae, Hippocastanoideae and Xanthoceroideae (Sapindaceae).

Ribosomal DNA internal transcribed spacers (ITS) and partial external transcribed spacers (ETSf) are popularly used to infer evolutionary hypotheses. However, there is generally little consideration given to the secondary structures of these small RNA molecules and their potential effects on sequence alignment and phylogenetic analyzes. Intergeneric relationships amongst three of the four major lineages in the Sapindaceae, the Dodonaeoideae, Hippcastanoideae and Xanthoceroideae were assessed by firstly, generating secondary structure predictions for ITS and partial ETSf sequences, and then these predictions were used to assist alignment of the sequences. Secondly, the alignment was analyzed using RNA specific models of sequence evolution that account for the variation in nucleotide evolution in the independent loops and covariating stems regions of the ribosomal spacers. These models and phylogeny drawn from these analyzes were compared with that from analyzes using 'traditional' 4-state models and previous plastid analyzes. These analyzes identified that paired-site models developed to deal specifically with stem structures in RNA encoding sequences more appropriately account for the evolutionary history of the sequences than traditional 4-state substitution models.